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ABSTRACT 


Recently-performed  analysis  for  herringbone  thrust 
bearings  has  been  Incorporated  Into  the  theory  of  the  vis¬ 
cous  screw  pump  for  Newtonian  fluids.  In  addition,  certain 
earlier  corrections  for  sidewall  and  channel  curvature  effects 
have  been  simplified.  The  result  Is  a  single,  refined  formu¬ 
la  for  the  prediction  of  the  pressure-flow  relation  for 
these  pumps. 

T*>  provide  a  self-contained  presentation,  some  of  the 
results  of  FIL  Interim  Report  I-A2049-31  are  also  given  here. 
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I.  INTRODUCTION 

Figure  1  shows  a  schematic  diagram  of  a  viscous  screw 
pump.  This  device  is  used  to  pump  relatively  small  quantities 
of  viscous  liquid  against  relatively  high  pressure  differences. 
It  Is  used,  for  example,  In  the  melt  extrusion  of  plastics. 

In  operation,  as  the  screw  rotates,  the  barrel  surface  drags 
by  viscous  shear  a  ribbon  of  liquid  between  the  •’flights**  of 
the  screw  from  a  feed  port  at  low  pressure  to  a  discharge  port 
at  high  pressure. 


Figure  1.  Schematic  Diagram  of  Pump 
(adapted  from  ref.  2) 

A  substantial  body  of  literature  pertaining  to  the  opera¬ 
tion  of  the  screw  pump  already  exists,  as  the  reader  who  pur¬ 
sues  tne  subject  will  discover.  The  references  cited  herein 
are  believed,  however,  to  give  a  state-of-the-art  picture  of 
the  theory  of  viscous  screw  pump  operation  with  Newtonian 
fluids.  The  purpose  of  the  present  work  has  been  both  to  gen¬ 
eralize  and  to  simplify  some  of  the  prior  theoretical  results. 
In  particular,  it  has  been  possible  to  Incorporate  almost 
directly  some  recent  analysis  performed  to  treat  herringbone 
thruct  bearings. 

The  simplest  theoretical  treatment  of  the  viscous  screw 
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pump  Involves  an  "unwrapping"  of  the  pumping  channel  to  create 

a  passage  between  parallel  planes.  Plow  over  the  flight  tips 

is  neglected.  The  channel  depth  Is  taken  as  very  small  compared 

with  Its  width,  so  that  Reynolds*  equation  can  be  employed  for 

computation  of  mass  fluxes.  A  simple  formula  relating  flow 

(5) 

and  operating  pressure  difference  can  thereby  be  obtained 

Successive  improvements  of  the  simplified  theory  are  pos¬ 
sible.  For  the  Isolated,  rectangular,  parallel-plate  channel 

the  effects  of  finite  width  ("sidewall  effects")  can  be  taken 
(4,6,2) 

Into  account  .  Also,  for  the  same  geometry  compatibility 

between  the  interior  transverse  pressure  ripple  and  uniform 

exterior  pressures  (at  inlet  and  discharge)  can  be  achieved 
(2,3) 

.  Without  unwrapping,  the  Isolated  rectangular  channel 

of  large  wldth-to-depth  aspect  ratio  can  be  analyzed  for  the 

(1) 

effect  of  channel  curvature  ,  If  sidewall  effects,  channel 

ourvature  and  entrance-exit  effects  are  neglected,  Reynolds* 

equation  can  be  solved  jointly  for  groove  and  ridge  film  thick- 

(3.7) 

nesses  .  Moreover,  variation  of  the  channel  depth  In  the 

(7) 

direction  of  flow  can  be  accommodated 

In  the  present  work,  a  Reynolds*  equation  treatment  for 
arbitrary  transverse  ridge  and  groove  shape  is  given  for  the 
pump  interior,  and  the  general  theory  for  the  corresponding 
edge  effects  Is  set  up.  The  edge-effect  analysis  is  numerical¬ 
ly  implemented  for  the  case  of  rectangular  grooves  and  ridges, 
an  exact  solution  being  obtained  in  the  case  of  an  Isolated 
groove.  Simplifications  of  the  treatment  of  sidewall  and 
curvature  influences  are  offered. 


2 


II  THE  INTERIOR  PRESSURE  SOLUTION 

The  starting-point  for  the  present  analysis  is  Reynolds’ 
equation  for  a  constant-property  Newtonian  fluid,  applied  to 


the  unwrapped  groove -ridge  combination  shown  in  Figure  2. 


Coordinates  and  Dimensions  Employed 


Thus : 


(2.01) 


Conversion  of  this  equation  to  the  skewed  coordinate  system 


(yi>x*)  gives j 


b  _  t,A  &  U  rr  f=  ~  f—  Xt  U  _  i» 
r  —  c*  J  A  Z  H  =  c 


Then: 


(2.03) 


?H3 


=  i  H'$ + H'#  -  2M£  &  if'  w  ( 

Asymptotically  (away  from  inlet  and  exit)  it  is  antici¬ 
pated  that  the  pressure  will  rise  linearly  with  distance 

along  the  groove.  Substitute  the  following  "interior"  solution 
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Into  eq.  (2,04). 


*  =  -t~~ —  +  fW  (2.05) 

A  /.  „  4 


An  ordinary  differential  equation  for  the  transverse  pressure 
ripple,  f(f),  is  thereby  obtained.  Thus: 


__  A  u3M  *  JH’ 

7f 

(2.06) 

One  integration  gives: 

H  —  H  H  +-  J  ;  jJ  a  comtont 

(2.0?) 

Or: 

«-3h" 

(2.07) 

Since  the  net  change  in  "f"  over  one  cycle  in  " f  " 

must  be 

zero,  "g"  must  satisfy  the  following  condition. 

O  -  ff*  4-  <*  -  (2.08) 

Therefore: 

W5  it  =  +  <2.°9) 

H3 /($)  =  -  ±  F($)  G($)  (2.io) 

Where :  ^ 

F(0=  2 S(Fhl JW- iP)dt  (2.11) 

0  O 

The  "interior"  pressure  solution  Is  therefore  given  by: 

kw+im/F  (2.12) 

For  rectangular  geometry  such  that: 

H  =  H, }  0^  £  *  r  :  H  =  HZj  Y<-  t  <1  (2.13) 
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the  pressure  solution  is* 


(£<»,/*)  7  '  9  J  Hm* 

Here  ^tC^)  1»  the  "sawtooth  function"  depicted  in  Fig.  3. 


Figure  3  Sketch  of  the  "Sawtooth  Function" 


III  THE  EDGE  CORRECTION 


Although  solution  (2.11)  is  an  exact  solution  of  the 
differential  equation  (2.04),  and  is  adequate  for  the  pump 
interior,  It  does  not  satisfy  conditions  of  pressure  uniformity 
at  f  =  0  and  =  1.  Therefore,  In  accordance  with  standard 
procedure  in  the  solution  of  linear  differential  equations, 
a  homogeneous  solution  is  sought  which  can  be  combined  with 
the  inhomogeneous  solution  to  satisfy  boundary  conditions. 

Let  cf(^if)  be  a  solution  of  the  following  homogeneous 
differential  equation: 

Boundary  conditions  on  "q"  appropriate  to  the  inlet  region 
are ; 

f(oj)  =  *  i  P(Q 


(3.02) 


f  V>,  f)  <  Constant,  as 
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Likewise,  boundary  conditions  appropriate  to  the  exit  region 


are  i 


f  (%}  t  ^  ^  }  d  ConStt+t)  f-X-^  °o- 


(3.03) 


Corresponding  functions  (f-  ?: )  (f- )  vanish  In  the 

pump  Interior,  but  possess  the  required  "ripple"  at  the  exit 
and/or  entrance.  Evaluation  of  specific  "  "  will  be  deferred. 


As  the  general  solution  satisfying  both  differential 
equation  and  boundary  conditions,  we  have,  then: 


The  constant  "ol"  In  eq.  (3.04)  can  be  found  by  taking  the 
difference  between  results  at  ^  =  0  and  £=  1.  Thus: 


*<>-*«»-  (ft;)  *  ?-  -  V-  ‘3-°s> 

In  view  of  the  definition  of  "IT"  and  the  fact  that  the  " 

are  Independent  of  (A/L),  It  can  be  seen  that,  as  (A/L) 
approaches  zero  with  finite  pressure  difference  between  en¬ 
trance  and  exlt,"ot"  Is  ascertainable  with  Increasingly  high 
precision.  In  other  words,  the  edge  phenomena  have  vanishing¬ 
ly  small  effect  on  the  general  linear  pressure  rise  within 
the  pump  as  the  groove-to-length  ratio  of  the  channels  becomes 
small. 

IV  FLOWRATES  ACCORDING  TO  REYNOLDS  EQUATION 

The  mass  flux  per  unit  transverse  distance  In  the 
x^-dlrectlon  Is  given,  according  to  Reynolds*  relations,  by: 
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(4.01) 


~£f&X 


Conversion  to  the  dimensionless  variables  f #  and  {  " 


gives : 


(4.02) 


Since  eq.  (2.05)  gives  an  accurate  evaluation  of  "IT"  in  the 
pump  interior,  and  since  the  total  throughput  must  be  invariable 
at  all  axial  locations,  this  equation  is  now  used  in  (4,02), 


sc  ui  cq.  (2.09)  then  yields: 

^  =  -  pc  Q+  H'fPtaSfi)  +  (fir*-  HPJ^(4.04) 


Therefore: 

- .  ///.  •  t-i3r  ..  ^ 

(4.03) 

Wc~p)  '  J 

Use  of  eq,  (2.09)  then  yields: 

">x 

The  total  mass  flow  per  groove-ridge  combination  is  found  by 
integrating  over  a  distance  Thus: 

Ai=  -JLsJ^LPlA  ^«0+  FT'FPtafp)*  (iP-HtroJ  (4.05) 

The  pressure-discharge  performance  of  the  screw  pump  can  now 
be  found  —  according  to  Reynolds'  equation  —  by  the  insertion 
in  eq.  (4,05)  of  the  pres sure -grad lent  value  of  "dC'  from 


eq.  (3.05). 

An  alternate  way  of  arriving  at  the  expression  for  "M" 
consists  of  calculating  the  flow  components  separately  across 
lines  AB  and  BC  shown  in  Fig.  4.  These  flows  combine  to  give 
"  M  " ,  of  course,  but  it  will  be  helpful  for  later  purposes  to 
possess  this  flow  breakdown. 
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Figure  4  Control  Region  ABC  Used  for  Secondary  Method 

of  Calculating  Mass  Flow 

The  mass  flux  In  the  Xj -direction  is,  according  to 
Reynolds*  approximations: 

^  '“-“i 

Components  of  the  mass  flux  vector  are  now  readily  found 
in  the  yf  and  y^ -directions  by  use  of  eqs.  (4.01)  and  (4,06). 

In  particular,  it  is  found  that: 

8  _ 

/ ft)  d/z  =-  PC  U&£*J3)  (x+H~*)A  (4  07) 

/  '  2  TF* 

and !  c  _ 

/ */•  =  f  -  lH  _<*  A  (4.08) 

8  *  * 

V  SIDEWALL  EFFECTS 

Viewed  in  the  ya -direction,  the  groove-ridge  cross-section 
within  the  pump  appears  like  that  in  the  sketch  of  Fig.  5 
Reynolds'  equation  cannot  take  account  of  the  fact  that  fluid 
velocities  must  vanish  on  the  sides  Pg  and  RS  of  this  configura¬ 
tion.  To  accommodate  such  boundary  conditions,  the  full  Stokes' 
equation  must  be  used. 

There  can  be  little  doubt  that  even  within  the  more  pre¬ 
cise  theory  of  Stokes' equation,  an  interior  region  of  the 
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Groove -Ridge  Cross-section 


y, 


pump  exists  where  velocity  patterns  f*re  established  and  where 
Mi  is  constant.  For  the  coordinate  system  (y,  ,y#  ,y? )  Stokes’ 
equation  for  becomes s 


Tyt >  3/3*  ^  V* 


(5.01) 


Although  this  is  a  simple  form  of  Poisson's  equation,  the  boun¬ 
dary  conditions  on  the  groove-ridge  corrugations  are  awkward. 

Since  usually  the  ridge  or  flight  film  thickness  is  much 
less  than  the  channel  or  groove  film  thickness,  examination 
of  the  solution  for  an  Isolated  groove  of  dimensions  {hff  x  s). 
Exact  solutions  have  been  given  by  a  number  of  authors  in 
several  different  forms.  Here  we  write: 


(jJ  ST 


with: 
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k*  _] 


*  -fr)  &  2  ' 

-|r  ftp-fa)  ^}J 


(5.  02) 


(5.03) 


The  series  terms  in  the  foregoing  equations  represent  devia¬ 
tions  from  the  results  of  Reynolds*  equation.  For  h  /s  ^3/4, 


q 


the  following  approximations  are  accurate  within  2.5%*  improv¬ 
ing  very  rapidly  as  the  depth-to-wldth  ratio  diminishes. 


Thus: 

z  V*  L 

7—  O.S42  8(j^~j  (5.04) 

*  ”  ^lyL  %  J 

1  -  0(>30z(2<!r)J  (5.05) 

A  derivation  for  the  flow  deficiency  along  the  sidewalls 
from  the  groove  bottom  to  some  arbitrary  position  " /3 "  is  given 
in  Appendix  A,  For  this  more  general  case,  expressions  (5.04) 
and  (5.05)  become: 

0."  i  KKs[l-  l5-061 

The  functions  "f j "  and  "f^"  are  shown  in  Figs.  6  and  7  respect¬ 
ively. 

The  foregoing  sidewall  analysis  is  strictly  applicable 
only  to  the  case  of  an  isolated  (or  closed)  rectangular  groove. 
However,  with  PQ/PT  taken  as  the  argument  (yj  /h<^  in  Figs.  6 
and  7#  the  flow  deficiencies  for  open  rectangular  grooves  are 
approximated.  This  particular  mode  of  approximation  will  lead 
to  somewhat  greater-than-actual  absolute  values  of  both  flow 
deficiencies,  since  some  of  the  flow  constraint  along  gT  is 
missing  in  the  case  of  the  open  groove. 
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VI  EFFECTS  OF  CHANNEL  CURVATURE 


It  is  usually  justifiable  to  treat  the  flows  in  the 
grooves  and  over  the  ridges  as  though  they  took  place  be¬ 
tween  parallel  flat  plates.  However,  it  is  possible  for  the 
radius  of  the  screw  root,  R^,  to  be  sufficiently  different 
from  the  radius  of  the  pump  barrel,  Rq,  so  that  the  unwrap¬ 
ping  process  can  introduce  some  error.  This  error  is,  of 

course,  principally  associated  with  the  groove. 

(1) 

Booy,  in  an  analysis  which  will  not  be  repeated  here, 
considered  the  effect  of  channel  curvature.  He  treated 
the  case  of  lnfinitesmal  film  thickness  over  the  flights, 
and  presumed  an  "infinite"  groove  (or  channel)  width. 

Using  Stokes'  equation,  Booy  showed  that  it  is  possible  to 
obtain  a  flow-field  solution  in  the  following  form. 

Radial  velocity  and  pressure  gradient  negligible. 

Axial  velocity  and  pressure  gradient  inter-related  as 
f  oIIovjs  : 


Tangential  velocity  and  pressure  gradient  inter-related 
as  follows: 


4T  =  A- jP,*L  tM. 
/-**(  r  /  2^.  0 -*’) 


(6.02) 


The  two  undetermined  pressure  derivatives  appearing  in 
the  foregoing  equations  were  determined  by  requiring  that: 

(a)  the  net  flow  normal  to  the  ridges  (or  flights)  be  zero. 

(b)  the  volumetric  throughput  be  appropriately  related  to 
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the  integrated  pressure  gradient  from  entrance  to  exit  at 
ny  fixed  transverse  position  in  the  groove.  (Note:  all  such 
integrated  gradients,  or  pressure  differences,  must  be  the 
same,  although  transverse  pressure  ripples  remain  to  be  can¬ 
celled  by  edge  corrections. ) 

In  the  case  of  helical  grooves,  the  assumption  of  "infin¬ 
ite"  groove  width  is  justified  only  when  the  groove  depth  is 
small  compared  with  the  screw  root  radius.  Accordingly,  it 


is  reasonable  to  restrict  consideration  to  those  cases  where 
the  ratio  (R^/Hq)  is  greater  than  0.5  In  Appendix  B  it  is 
proved  that  a  virtually  "flat-plate"  formula  results  if  a 
"mean"  radius  is  selected  two-thirds  of  the  way  from  the  jt 

to  the  barrel  radius.  The  helix  angle  at  this  mean  radius 
is  taken  as  . "  The  simplified  formula  is: 


a  =  s 


il±jkl  l~  14.4 _  Jk.  2fc)l-  *--hs 

+  #>*»  L2-  ^  yj'  =  **. 


(6.03) 


Tills  formula  shows  that  within  the  limits  designated,  the  ef¬ 
fective  groove  width  varies  by  about  plus  or  minus  8#.  Fur- 

shows  that  the  curvature  correction  virtually  vanishes 
for  a  groove  angle  of  45  degrees. 

Booy\s  solution  of  Stokes*  equation  does  not  yield  zero 
local  velocity  normal  to  the  groove  edges  -  only  zero  net 
flow  normal  thereto.  This  same  tyre  of  average  conformity 
to  the  velocity  boundary  condition  would  be  provided  by  a 
solution  of  Reynolds’  equation  in  the  case  of  a  more  shallow 
groove.  Hence  sidewall  corrections  are  in  order.  Rigorously, 


such  corrections  depend  on  the  radius  ratio, o(.  However, 


If  corrections  on  corrections  are  dismissed  by  reason  of 
their  second -order  smallness,  the  analysis  of  Section  V 
Is  applicable. 


VII  fc DGE  SOLUTIONS  FOR  RECTANGULAR  GROOVING 


At  this  point  we  revert  to  further  consideration  of  the 
entrancc-exlt  edge  problem.  It  suffices  to  treat  this  aspect 
of  the  pump  problem  with  Reynolds'  equation.  Moreover,  atten¬ 
tion  will  be  directed  to  the  special,  but  extremely  important, 
case  of  rectangular  grooving  with  equal  groove  and  ridge 
width. 

The  analysis  of  Section  III  demonstrates  that  the  princi¬ 
pal  result  required  from  solution  of  an  edge  problem  is  the 
value  of  q^ .  Some  Information  concerning  this  value  is  pro¬ 
vided  by  the  mass-content  rule.  Thus,  integrate  eq,  (3.01) 
over  one  groove-ridge  cycle; 

(i f&fW  -  fTpr-° 

°  o'  o  (7.01) 

Integration  of  the  second  term  by  parts  results  in  some  cancel¬ 
lation,  and  the  following  simpler  equation. 

/  1 

/t)£>  /A/Vf  +  /r ft = °  (?-o2) 

O  Q  * 

This  equation  can  be  integrated  once  with  respect  to  "if," 

The  integration  constant  must  be  zero  if  uniformity  is  to  be 
achieved  at  infinity.  Therefore; 


(7.03) 
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In  the  care  of  rectangular  grooving,  the  fllm-thlckness 
derivative  becomes  a  delta-function  at  /*  =  0  and  J  =  Y  . 

Thus :  | 

o  (7.04) 

It  would  appear  that  further  progress  cannot  be  made  along 
there  liner  for  the  care  of  arbitrary  groove  angle,  since 
q*Q  and  depend  on  both  and  the  f llm-thlcknerr  ratio, 
^1^2*  pPRcia^  care  of  ^  ~  ir/2,  It  Is  seen  from 

eq.  (7.03)  that  In  all  cares: 

/ 

Jtl’f  —  Constant  (7.05) 

o 

With  rectangular  grooving  (and  any p  )  q'  must  become  the  saw¬ 
tooth  function  along  £ =  0.  (See  eq .  2.14  and  Fig.  3)  For 
convenience,  assign  unit  amplitude  (rather  than  K* )  to  this 
function.  Then  it  is  easy  to  shov;  from  eq.  (7.05)  that: 

=  0.5 

•© 

For  arbitrary  "g"  a  solution  for  "cf  "  probably  cannot  be 
obtained  In  closed  form,  even  for  rectangular  groover  and 
ridges.  However,  for  an  Isolated  rectangular  groove  (zero 
ridge  or  flight  film  thickness),  the  author  has  been  able  to 
obtain  such  a  solution.  Details  of  the  analysis  are  given  in 
Appendix  C.  The  procedure  involves  transforming  by  conformal 
mapping  problem  A  of  Fig.  8  into  problem  B.  With  constant  film 
thickness,  "q”’  satisfies  Laplace* r  equation  In  the  ( /,  tyz) 
plane.  By  virtue  of  the  properties  of  conformal  mapping,  it 
will  also,  then,  satisfy  Laplace’s  equation  in  the  (U  ,&*) 
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Figure  0  Conformal  Transformation 
plane.  Zero  normal  derivative  along  the  groove  sides  Is  also 
preserved . 

Now  the  average  value  of  at  any  level  v  =  const.  Is 
also  "q£"  In  the  (u,v)  plane.  However,  both  planes  possess 
the  same  "q^.  ’*  Therefore,  the  determination  of  "q^"  Is  reduced 
to  an  Integration  along  v  =  0.  The  final  result  Is: 


|£=  fad)  -  ^(r)}  (7.06) 

where  ip  the  Digamma  function,  tabulated  (e.g.)  in  ref,  7. 

In  the  case  of  rectangular  grooving,  the  asymptotic  value 
of  "Ym  at  can  be  found  from  that  at  the  inlet. 

The  Inlet  and  exit  conf lgurat Ions,  as  shown  In  Fig.  9,are  help¬ 
ful  In  understand lng  the  relation.  The  pressure  ripples  which 
serve  as  boundary  conditions  for  the  respective  "q^'s"  are 
shown.  Now  the  "response"  to  the  pressure  ripple  of  the  Inlet 
conf  lguratlon  Is  "cf4 , "  The  response  to  a  uniform  pressure  of 

Oo 

unity  would  be,  of  course,  just  1.0.  Because  of  the  linearity 
and  homogeneity  of  the  differential  equation  for  "q" ,  the  res¬ 
ponse  to  (1.0  minus  the  inlet  pressure  ripple)  would  be  1  -  tf4 
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Figure  9  Sketch  to  Assist  Understanding  of  Inter-relation 

of  Inlet  and  Exit  '  s 

However,  as  can  be  seen  by  examining  the  ripple  for  the  exit 
configuration,  this  newly-created  compound  ripple  at  the  inlet 
is  precisely  the  same.  Therefore: 


(7.0?) 


In  particular,  for  the  Isolated  rectangular  groove: 


('-£)-  ^(i)}  (7.08) 


Reference  to  eq.  (3.05)  that  the  difference  r- 

is  the  quantity  needed  for  over-all  pump-performance  calcula¬ 
tions.  Accordingly,  this  is  the  quantity  presented  in  Fig.  10 
as  a  function  of  groove  angle  and  f ilra-thickness  ratio.  The 
groove  and  ridge  widths  are  taken  as  equal.  The  curve  for 
h2/hi  =  0  corresponds  to  eq.  (7. 06).  The  curves  for  other 

film-thickness  ratios  were  obtained  by  use  of  the  finite-element 

(B) 

numerical  approach  developed  at  The  Franklin  Institute  .  For 
this  data  the  author  is  much  indebted  to  Dr.  T.  Y.  Chu. 

Appendix  D  gives  the  numerical  values  on  which  the  curves  of 
Flc.  If-  i.;r«  bused.  Cross-plotting  was  employed  with  the 
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Figure  10  Pressure  Correction  Fade 

Rectangular  Groc 
(Versus  Groove  Angl 
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(1.0) 


X  —  MUVDLRMANH'S 
Aa/ALOQ  results 


ANGLE ,  (3 


**  '  -vvi 


/or  Effect  with 


assumption  of  linearity  in 

The  asymptotic  levels  of  "  can  be  found  from  Fig.  10, 
but  no  Indication  is  provided  concerning  the  distance  along 
the  groove  or  channel  required  to  attain  them.  A  reasonable 
estimate  can  be  formed  by  further  consideration  of  the  case  of 
the  Isolated  groove  with  In  this  particular  case,  a 

standard  Fourier  treatment  yields: 


Ah  C*4,(nvu) 

hs.Q 


-  MTV 

e 


(7.09) 


Here  AQ  =  y  =  .  Deviations  from  uniformity  approach  zero 

exponentially,  with  the  least  rapid  attenuation  associated  with 
the  first  mode,  containing  the  factor  €  .  When  v  =  1, 

this  factor  has  already  been  reduced  to  only  4.3$  of  its  value 
on  the  boundary  v  =  0.  Translated  into  a  more  general  con¬ 
clusion,  this  result  means  that  end-effects  persist  for  only 
about  one  channel-width  into  the  channel.  In  the  case  of 
oblique  channels,  measurement  towards  the  interior  should  be 
made  from  the  apex  of  the  obtuse  angle.  To  achieve  asymptotic 
conformity  within  the  pumn  interior,  there  is  a  minimum  axial 
length  required.  It  is  given  by:  (approx.) 


-=  >:  ZA**/8  + 


(7.10) 


VIII  COMPARISONS  WITH  OTHER  INVESTIGATIONS 

In  order  to  develop  a  suitable  pump-performance  formula, 
it  has  been  necessary  to  repeat  some  material  presented  by 
earlier  investigators .  Here  an  effort  will  made  to  indicate 
the  aspects  of  novelty  of  the  present  work,  to  inter-compare 
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with  earlier  results,  and  to  report  various  checks  for  accuracy 
which  were  made. 

Accurate  account  for  cross-flight  leakage  (within  Reynolds' 

(9)  (10) 

lubrication  theory)  was  made  by  Whipple,  Vohr  and  Pan  and  by 

others  treating  herringbone-Dattern  bearings,  as  vrell  as  by 

(3) 

Mul jdermann  in  work  cited  earlier.  However,  the  derivation 
given  in  flection  II  is  believed  to  be  new,  and  is  simpler  than 
the  matching  of  solutions  at  groove-ridge  Interfacial  discontin¬ 
uities.  For  the  special  case  of  rectangular  grooving,  eq.  (2.14) 
gives  Mui jdermann ' s  "approximate  pressure  distribution." 

Lxact  solutions  for  established  laminar  flows  have  a  ven¬ 
erable  past.  "Pressure-f low"  expression  (5.03)  for  a  rectangu- 

(4) 

lar  flow  cross-sect  ion  was  given  by  Carley  and  Strub  as  a 
form  preferable  to  that  given  by  Boussinesq  —  although  mathema¬ 
tically  equivalent.  In  the  same  way,  expression  (5.02)  for 
"drag  flew"  is  somewhat  to  be  preferred  to  equivalent  series 
given  in  refs.  (1),  (4),  and  elsewhere.  In  each  case,  rapid 
series  convergence  with  attendant  computational  simplicity  is 
obtained  by  treating  the  sidewall  flow  deficiency  as  a  residual, 


boundary-layer  effect,  brought  out  further  by  the  approximate 
expressions  (5.04)  and  (5.05).  The  partial  deficiency  functions 
in  eqs.  (5.06)  and  (5.07)  capitalize  on  this  simplicity. 

Booy * s  analytical  results  for  the  influence  of  channel 
curvature  are  reproduced  by  eq  (6.03)  in  the  (hG/R0)  range  of 
most  interest.  With  virtually  the  same  degree  of  approximation, 


his  flow  factors  are  given  by:  (his  notation) 


i 

i 

i 
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(B. 01 ) 


Values  from  eqs.  (8.01)  agree  adequately  with  those  from  Booy's 


Figs.  3  and  4 

of  ref.  (1). 

The  very 

worst  agreement,  corres- 

ponding  to  a 

radius  ratio 

of  2. 0  and 

extremes  of 

groove  angle 

is  shown  below. 

F  (8.01) 

DC 

F  (Booy) 
DC 

F  (8.01) 
PC 

F  (Booy) 

PC 

ii 

1.08 

1.07 

1.242 

1.308 

=  60 

0.722 

0.693 

0.833 

0.84? 

Agreement  improves  rapidly 

(3) 

with  diminishing  "£  " 
(2) 

• 

Both  Muljdermann  and  Booy  have  dealt  with  the  edge 
effect  at  the  inlet  and  exit  of  an  isolated  groove.  Muij- 
dermann's  electrical-analog  results  are  shown  on  Fig,  10. 

His  results  agree  with  the  present  exact  results  within  20 %9 
and  the  agreement  is  considerably  inside  this  limit  for  most 
measurements.  It  should  be  recalled  that  any  deviations 
here  are  in  a  correction  term.  A  less  demanding  comparison 
Is  possible  with  Booy's  finite-difference  computations  of 
his  pressure-flow  correction  factor,  Fp^ ,  This  factor  satis¬ 
fies  the  following  relation:  (Booy's  notation) 


£-  t- 


re 

Values  from  this  equation  are  imperceptibly  different  from 
those  presented  in  Booy's  Fig.  12,  where  curves  for  ft  -values 
of  arctan(l/3il/2,l)  were  presented. 

For  the  case  of  rectangular  grooving  with  H2/H^  4  0, 
there  have  hitherto  been  no  accurate  corrections  available. 
Mui Jdermann  tentatively  offers  a  formula  based  on  heuristic 
reasoning.  He  expresses  all  pressure  corrections  in  terms  of 


(8.02) 
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the  pressure  correction  for  an  isolated  groove.  Thus: 


Although  this  formula  gives  correct  limiting  values  at 
H?/^l  =  0  an(J  =  It  It  would  appear  to  be  considerably 

in  error  for  some  intermediate  values.  For  example,  with 
/i  =  15'  and  (H,,/^)3  =  0.06,  Fig.  10  gives  Aj^  =  0.43, 
whereas  eq .  ( 8 . 03 )  gives : 

1  -  0.06 

0.79  -  =  0.70 

1  +  0.06 

Actually,  the  curves  in  Fig.  10  for  ^2^1  ^  0  must  be 
regarded  as  somewhat  tentative,  inasmuch  as  the  finite-element 
calculations  were  not  "converged"  by  running  wltn  diminishing 
grid  size.  Checks  with  the  exact  conformal-mapping  solution 
warrant  an  error  estimate  of  less  than  3 %• 

IX  PUMP  PERFORMANCE.  FORMULA 

The  individual  results  of  the  preceding  sections  will  now 
be  combined  to  yield  the  overall  relation  between  pressure  dif¬ 
ference  and  quantity  pumped.  The  development  will  be  for  one 
groove-ridge  (channel-flight)  pair,  and  it  will  be  assumed 
that  the  channel  width  has  been  adjusted  for  curvature  prior 
to  "unwrapping"  by  use  of  the  factor  in  eq.  (6.03). 

With  reference  to  Fig.  4,  the  flow  across  AC  is  equal  to 
the  flow  across  AB  plus  the  flow  across  BC.  Low  Reynolds' 
equation  is  deemed  satisfactory  for  the  computation  of  all 
flows  excent  that  alon/r  the  groove.  (The  transverse  flow  over 
AB  Is  small,  in  any  event,  and  the  flow  along  the  ridge  in- 
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volves  a  very  small  film  thickness, )  The  sidewall  flow  defic¬ 
iency  in  the  groove  is  given  by  the  second  (subtractive)  terms 
in  eqs,  (5.06)  and  (5.07).  Thus  the  total  flow,  including 


sidewall  effects,  is  given  by: 


Q  =  Q  -  Q  (9. 

Rey  def 

The  flows  will,  of  course,  be  evaluated  in  the  pump  Interior. 


(9.01) 


Again: 


&=  Q_  -j^sF  %s  &)7+l 

Ttr  L  z  rJ 


(9.02) 


The  factors  F^  and  Fp  will  ordinarily  be  given  well  enough  by: 


(9.03) 


The  required  pressure  derivative  can  be  put  in  dimensionless 


form.  Thus: 


-  ft} 


19.04) 


When  this  derivative  Is  used  on  the  Interior  pressure  solution, 
given  by  eq.  (2.14),  and  Inserted  Into  eq.  (9.02),  the  result 


Combination  with  the  Reynolds'  flow  now  gives: 


-  cjOjCaCf 
2/ 


(9.05) 


t  /^=  +  Hs  -  yfrt/0  Hi  Fpj  J  ( 9 . 

+  -  H  +  r^> 
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Further  straightforward  algebraic  rearrangement  gives: 

Let  us  turn  now  to  the  pressure-d 1 f ference  expression 
given  by  eq .  (3.0$).  The  "qf^s"  of  this  equation  must  corres¬ 
pond  to  the  ripple  term  In  eq .  (2.14),  since  rectangular 
grooving  Is  nov  assumed.  For  unit  sawtooth  ripple,  Fig.  10 
give-?:  —  l  5e  Zi/e,*-*. 

Hence,  for  rectangular  grooving,  eq .  (3.0.5)  becomes: 


0?) 


TO)-  Tr(»)  = 


(A  a*.p) 

where  the  amplitude  is: 


(dmf>L/^uJie )  Aj?c 


(9.08) 


(9.09) 


AmpHtndt  ~  +  «  (H?~  H-0j 

Again,  some  algebraic  rearrangement  gives: 

cM/>  =  «  -  4^  f c((fP-  u;3)  +  4-Wj  (9 

6*1,1/,  Hup  0  L  H'*  J  ' 

hlirninatlon  of  the  pressure-grad  lent  term  "oc."  between  eqs. 
(9.10)  and  (9.0?)  yields  the  following  final  pump  performance 


10) 


AV 


Here  : 


Thus:  (with  c  =  hG , 

so  that  E1q  =  H-J  =1) 

=  / 1  —  A^r  Stl  I  ±  S*.  —  & 

l  L  rr~r  \ 

p9.11) 

AP=  cAp 

.  Q 

(9.12) 

La  Uj  (CLa^ 

■?”  i h  ; 

Sz  =  !-  ir- 
FF 

(9.13) 

=  =  +  W’U/i  -  yfp  tJifi 

(9.14) 
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A  NUMERICAL  EXAMPLE 


Figure  11  Dimensions  of  Experimental  Screw  Pump 

The  dimensions  of  an  experimental  screw  pump  tested  at 
Columbia  University  are  shown  above  In  Fig.  11.  Data  for 
this  test  were  kindly  supplied  the  author  by  Prof.  D.  D.  Fuller. 

The  pump  possesses  two  single-thread  pumping  sections  working 
In  opposition.  The  channel-grooving  is  in  the  barrel.  Per¬ 
tinent  characteristics  are  listed  below. 

Pitch  of  Threads:  four  per  inch 
Shaft  Diameter:  2,1192  Ins, 

Barrel  Diameter:  2.123B  ins. 

Width  of  Channel:  0.125  Ins. 

Width  of  Flights:  0.125  ins. 

Depth  of  Channel  Grooving:  0.050  Ins. 

The  corresponding  film  thicknesses  are: 

Ridge  Film  Thickness:  hR  =  radial  clearance  =  0.0023  ins. 

Groove  Film  Thickness:  hn  -  0.0023  +  0.050  =  0.0523ins. 

G 

Also:  Hg  =  Hx  =  1  and  H2  0. 0023/° .  0523  =  0.044 

Unwrapping  of  Pump  Channel 

Although,  in  this  instance,  the  effect  of  channel  curvature 
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is  very  small,  for  the  sake  of  Illustration  Its  magnitude  will 
be  comnuted. 

=  0.5  *  2.1192  +  (2/3)*  0.0523  =  1.09^5  ins. 

Therefore : 


0.250 
2ir*  1.09' 


^  0.0364  =  ta 


In  addition: 


e  =  h,./n  =  0.0523/(0.5*2.1192  +  0.0523)  =  0.047 

Th?  effective  groove  vldth  Is  now,  according  to  eq.  (6. 03), 
sRff  =  0.125  (1  +  0.047/6)  =  0.126 
Sidewall  Factors: 

For  determination  of  the  factors  F^  and  Fpf  we  have: 
hG/s  =  0.0523/0.126  =  0.415 
and:  y3/hG  =  O.O5/O.0523  =  0.956 

V'ith  the  definitions  in  eqs.  (9.03)  and  the  use  of  Figs.  6 


and  7 1  >’o  obtain: 


Fd  =  0.415  *  0.516  =  0.214 
Fp  =  0.415  *  0.624  =  0.259 
Film-Thickness  Constants : 

The  following  different  averages  of  film  thickness  are 


required . 

H“  =  0.5  (1  +  0.044)  =  0.522 
i?=  0.5  +  0.5*(0. 044)3  =  0.500 
H*=  0.5  +  0.5/(0.044)2  =  260 
?*=  0.5  +  0.5/(0.044)-^  =  5365 


Evaluation  of  the  S^: 

The  four  S-constants  can  now  be  evaluated.  Thus: 
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Sj  =  1  -  1/5365  =  1.000 

52  =  1  -  260/5865  =  0.955 

53  =  .522  -  260/5865  -  .5*. 214  =  0.370 

=  1/5365  +  ,5*(.0364)2  -  ,5*.259*(  .0364)2 
=  6.62*10"^ 


Final  rump  Formula : 

Direct  sul't. itu! \>>n  Vt:»  *q.  (9.11)  nov;  gives  the  final 
formula  relating  operating  pressure  difference  and  quantity 
of  liquid  being  pumped.  (Apcorr  Is  read  from  Fig.  10.)  Thus: 

AP*  +  *21 5  =  h  -  .5*x.252*,g75*i*il  f.,2 20  .  A  +  *251  -  _  —9HL—X 

"  1  l  1  JIO25T0-*  1  6,62*10"*  *1/ 

Or: 

AP*  =  488  -  1320  Q*  (10.01) 


It  will  be  observed  that  the  numerical  constants  In  the  above 
formula  depend  on  the  geometry  only.  The  speed  and  the  fluid 
viscosity  appear  solely  In  the  dimensionless  variables  and 

Q* ,  and  are  assignable  over  wide  ranges. 

Comparison  with  Experiment : 

In  an  experiment  conducted  with  the  pump  of  Fig,  11,  the 
following  measurements  were  made: 

Shaft  Speed:  3000  rpm 
Shut-off  Pressure:  *+2.5  psi 
Fluid  Viscosity:  4. 57*10**  lb-sec/in2 
The  viscosity  value  was  Inferred  from  a  viscosity-temperature 
curve  for  the  oil  with  measured  temperature  used  as  argument. 

A  prediction  of  the  shut-off  pressure  is  simply  made  with 
the  use  of  eq.  (10,01)  and  the  definition  of  Ap*  In  eq.  (9,12). 


28 


The  velocity  at  the  outer  radius  of  the  groove  (Kith  the  shaft 
treated  as  rtatlonary)  lrs 

U1  =  (3000/60) *3. 142*2*1. 112  =  348  ln/sec 
Then,  since: 

Af>  »  L  (A  £p*  (10.02) 

c 

the  prediction  of  shut-off  pressure  is: 

4,  g  488  -  56.3  psl 


In  a  second  experiment  with  the  same  pump,  the  data  were 
as  follov\s: 

fhaft  Speed:  3000  rpm 
Flow;  0.525  inVsec 
Pressure  Difference:  31.5  psi 
Fluid  Viscosity:  4,57  lb-sec/ln^ 

With  the  given  valuer  of  pressure  difference  and  viscosity; 

ZiP*  =  248  and  Q*  =  (488-248 J/1320  =  0.182 
Reversion  of  the  definition  of  Q*  in  eq.  (9.12)  gives; 


tteS-Lct  cT=  UcA-  <3* 

*  1  Z 


(10.03) 


Therefore : 


Q  = 

Compared 
O.525  In3/sec 
tlons  of  56.3 


( 3^8/2)*.  0523* ( .250/1)*, 182  =  0.415  inVsec 
(for  two  Dumping  sections)  ^B^O^n^seg 
with  the  measured  values  of  42.5  psl  an d 
for  the  shut-off  pressure  and  flow,  the  predic- 
psi  and  O.83O  ln-Vsec  are  optimistic.  However, 
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If  only  cross-flight  leakage  had  been  allowed  for,  with  no 
consideration  of  channel  curvature,  sidewall  effects,  and  Inlet 
and  exit  corrections,  the  values  would  have  been  62.8  psi  and 
1.12  in^/sec.  Clearly,  conslderat Ion  of  these  additional  ef¬ 
fects  greatly  improves  accord  between  theory  and  experiment. 

An  examination  of  possible  sources  of  the  residual  dis¬ 
crepancy  leads  to  the  conclusion  that  the  viscosity  values 
employed  in  the  calculations  are  too  high.  Such  values  would 
arise  through  measurement  of  temperatures  lower  than  actually 
existed  in  the  oil  films  —  an  eventuality  which  Frof.  Fuller 
suggests  as  highly  likely  in  the  foregoing  experiments.  In 
fact,  it  is  notevjorthy  that  if  the  viscosity  valuer  used  in 
the  theoretical  predictions  are  arbitrarily  lowered  by  2^1 , 
the  predicted  shut-off  pressure  becomes  the  same  as  measured, 
and  the  predicted  flow  is  a  mere  0.025  ln3/sec  higher  than 
measured . 
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APPENDIX  A 


SIDEWALL  ErrECTa 


It  Is  convenient  to  solve  eq.  (5.01)  In  terms  of  dimen¬ 
sionless  variables: 

Y* 

Then  : 


4' 


=  -A 


T'l  ~  ^  2 


* 


(A. 01) 


A  particular  solution  Ip  provided  by  heynoldr'  equation.  Thus: 

(A. 02) 


v’i 7-7) 


A  homogeneous  solution  for  "v^"  must  now  be  found  which  vanishes 
on  ^7=0  and  1.  This  solution  is  formed  from  the  series: 


/^N  = 


Z_ 

**/ 


where : 


f\(-r)  =  C0tl{*r(r-*fc)} 


«*cg) 


On  T*  =  0  and  1 : 

A3 


Then : 


>m.  (XTTJ)  ,  (f-  j) 


C*  H  +  l 


-  $  %  £  ^ 


(A. 03) 


(A. 04) 


(A. 05) 


(A. 06) 


Since  vH  muet  be  rubtracted  from  vp  to  give  the  true  Vg,  It 
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represents  the  velocity  deficiency  occasioned  by  the  sidewalls. 

The  flow  deficiency  per  unit  distance  along  one  sidewall 

is  given  by: 

_X. 

O  K  +1 


bf'V*  ~  ^ 

O  Ac  */ 

'4‘t  &  2{l=w^0av 


(A. 07) 


The  total  flow  deficiency  for  one  sidewall  Isj 

V/  /* 


o  o 


*-  •• 


(A. 08) 


When  this  last  result  is  multiplied  by  two ,  and  subtracted 
from  the  usual  Reynolds  flow,  eqs.  (5.02)  and  (5.03)  are  the 
consequences . 

As  stated  in  the  main  text,  when  h q/s  is  less  than  3A» 
sidewall  effects  are  essentially  confined  to  the  channel  ends. 
For  this  case,  all  the  hyperbolic  tangents  in  the  foregoing 
series  can  be  taken  as  unity.  Accordingly,  let  us  make  this 


simplification  and  compute  next  the  partial  flow  deficiency: 


:a. 09) 


Then: 


W  AT  I 

*Z9t>- 


(A. 10) 
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and  : 


£(%)-  g'JTOztoyp- (a.  11 ) 

The  series  for  f2  converter  exceedingly  rapidly,  and  can 
be  used  for  computation  directly,  with  the  results  shown  in 
Fig.  7.  The  series  for  f-j  is  best  converted,  as  follows. 

From  ref.  11,  n.  139.  we  have: 


Oo 


JL  Mt.  (>—  Zb  c**x  +>•  *)  =  x  a* 

*=/ 


(*x) 


(A, 12) 


retting  r  -  -1,  and  performing  two  integrations,  one  finds: 

x  x 

J  J &,{Z  C*£)}  dx  dx  =r  <WCW  (A.  13) 

0  o 


k=t 


Again,  from  ref.  11,  p.  123: 

zx-t 


J&H,(cO€x)  =  —  J/*  Z  (z  —  l)  “B*  X 

*  Oik)! 


(A. 14) 


k»i 


where  the  are  the  Bernoulli  numbers,  a  fev;  of  which  are 
given  below. 


B-,  =  1/6;  B2  =  1/30;  B3  =  1/42;  B4  =  1/30;  B5  =  5/66  . 

The  Integrations  of  eq ,  (A.13)  are  readily  performed  with  the 

new  series,  and  the  result  is: 

mo 

Zc-O^fj- *«**)}=  (A-15) 

k *  Z  ZK  (Z*+3.)j 

Finally : 


jL  ?l_ 

S7Lo1  IbIZfo 


Thin  function  Is  rraphed  In  Pip,  6. 


(A. 16) 
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APFENDIX  B 


CHANNEL  CURVATURE  EFFECTS 

(1) 

Booy  gives  the  following  formula  for  the  flow,  Q,  In 
term?  of  the  overall  pressure  difference,  Pq,  the  radius 
ratio,  o(  =  R^/Rq,  and  other  variablesf  Thus: 

Q  =  tz a/ S(i«) -  T(4„«)3L  < b. oi ) 

Here : 


f.L 


J=  Z*F(«)K(«)t^A  .  T—  nrFMCtitJy. 
F(<)-ZQC*)  FW-2GM 


( B. 02 ) 


m,~  xO 


(B.03) 


(B.04) 


!<(*)=  1  + 

0-*V 

In  succeeding  manipulations ,  use  will  be  made  of  the 
fact  that: 

h  —  ^0— 7^  ;  ~R  ielu,<f>  sr  ~W0  icu*  <j>o  for  any  R 


(B.05) 


(B.06) 


^Booy's  notation  used  in  this  Appendix.  A  factor  of  4  has  been 
introduced  Into  the  denominator  of  the  second  term  of  (B.01) 
to  correct  for  misprint  in  (1). 
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Further,  y2  is  chosen  as  distance  parallel  to  the  flights  at 
radius ,  R,  and  V 2  a?  the  corresponding  component  of  barrel 
velocity , 


Algebraic  rearrangements  and  substitutions  give  succes 


sively : 


<2_  \zrFC*)'ti^4ia 


zjrKNkM  _  i  3_ 

z  (V3)  tyi-l 


( B. 07) 


Q—  rrKKliyJ  F(*)kb) 


z*F.  N((-j)  Tl 

z  (Vi)  KW 


(B. 08) 


Q  —  .zfrff  &*fi/uc4  Fbdkc*) 

0-4){m-zmhi+m }  L 


vJj  _  fzel  6-*)} -L, 
Z  l  (%)  Kfrj  I 1 


(B.09) 


(%)  Kfrj 

Equation  ( R.  0°)  1  valid  for  any  "R."  Let  us  choose 
H  to  effect,  as  much  as  possible,  correspondence  with  standard 
flat-plate  formulas.  Thus,  take: 


Tm  V*&(«)  0-«)  m  Yo  (,-d)\  b3  (  b.  io ) 

(V3 )  KC«) 

Or: 

&)6-*FK60 

~  - GOO -  (E-n> 


Because,  lr.  mort  practical  circumstances ,  ’V"  can  never 
deviate  much  from  unity,  It  Is  helpful  to  vrlte:  tar  t-vL 
and  to  expand  the  various  functions  of  "«t  '  In  terms  of  " t  ", 
Thu  - ; 


-4^(i_±)  .  £$()=  */(nl)  .  K(«)=  e(i+f) 


(B.12) 
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(B.13) 


Formula?  (B.12)  are  very  satisfactory  approximation?  for 
the  practical  range  of  V" .  For  example,  even  when  hj /H^  =  1/2 
( =  1/2),  th .  exact  and  approximate  valuer  of  these  functions 
differ  little.  (Tee  the  table  below.) 

Fxact  Approx. 

F(.5)  -.126  -.125 

C(.5)  0 . 108  0.104 

K( . 5)  0. 533  0.543 


Ure  of  eq.  (B.ll)  in  eq.  (B.09)  rives: 

When  the  € -approximations  are  used,  the  following  simple 
equivalent  flat-plate  formula  results: 


(3.14) 


a 


_  {  jr  £/*.  !(/+{■) 


(B.15) 
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APPENDIX  C 


SINGLE-GROOVE  EDGE  CORRECTION 

The  mapping  function 

*  6  . 

r  £-'/  7^  *(*  / 

w  =  Jz  (*-•)  e  di  =  f  +  iy 


converts  the  real  axis  of  the  e-plane  into  the  form  shown 
below  in  the  w-plune 


x+  i  y  w  -  (i-{)+  ‘  ^ 


When  (3  =  —  we  take  f  S'  ^  }  ^  ^  tt 

Along  y  =  c 


Now  in  the  1/-)  plane 
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'  I 

-  _  jXO- 

=  i  - i* 


[Note 


=  1- 


:  f>r) 


Next  Let 


x-_il  •  t-JX 

|+  ia  J  'Vi-x 


then: 


/V  * 

?-  =  J 


te'(t)  t  *Ji 

t-t-i3- 


Write  tan  (t)  in  the  following  integral  form. 


Thus : 


1  -B 


••  ?  -  1-  2i±tf'>{z  {—L  dt  h 

/-  IT*  J  }(HiV>(n-t‘> 

J!r*0  t^O 
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Resolution  into  partial  fractions  gives: 

£- 1  -  tyj'iyzsh)  -i4  sW 1 


¥ 

dids 


Use  is  made  of  the  infinite  integral 


I 


i  jit  _  ch  r» 

cn-tr"  r(mto 


and  the  reflection  formula  for  the  Gamrru:  function: 


Then : 


n\ 

2,  1  _  -2  f  itl  1 1  ds  \  m  =  I  - 


The  above  definite  integral  converges,  but  not  with 


m 


the  terms  1,  5  treated  separately.  Accordingly,  the  behav¬ 
ior  of  an  integral  with  the  denominator  "shaded"  to 

/  M,:r 

(J —  S  )  will  first  be  examined. 


i-  =  i 

*  /.  ,X\'— 

__  — 

r0£') 

PC*) 

z 

r(%'* 
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Let  0  =r  h+J  ,  #  and  examine  "I"  as  6  -•  o. 
2* 


I-  ± 


Hence : 


e  P(e)  [i 

p(i)  r(«+ "t'j  -  r  ("40  rce*i  >] 

Ttfri)rfe+'®). 

- 

© 

J 

1 

rrft)pr^')_  rps-’jpfeil 

rd)  rC^1) 

L 

J 

i  _ 

1 

■’T'-l)  _  p'ft>~ 

T 

] 

’(»-£)  rm 

-  ‘fee  {«l>0-£)  - 


’More  conveniently  for  computation: 


f-  -  1  -  V  { iK*_£)  -  ^  ] 
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